Phase diagram of the SO(n) bilinear-biquadratic chain from many-body entanglement 
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Here we investigate the phase diagram of the SO(n) bilinear-biquadratic quantum spin chain by 
studying the global quantum correlations of the ground state. We consider the cases of n = 3, 4 
and 5 and focus on the geometric entanglement in the thermodynamic limit. Apart from capturing 
all the known phase transitions for these cases, our analysis shows a number of novel distinctive 
behaviors in the phase diagrams which we conjecture to be general and valid for arbitrary n. In 
particular, we provide an intuitive argument in favor of an infinite entanglement length in the 
system at a purely-biquadratic point. Our results are also compared to other methods, such as 
fidelity diagrams. 
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I. INTRODUCTION 

The recent ability to experimentally manipulate ultra- 
cold atoms in optical lattices has opened the possibility to 
simulate a number of strongly-correlated systems In 
fact, it is now conceivable to recreate the physics of some 
quantum spin systems of relevance in condensed matter 
physics . This great advance in quantum simulation has 
opened the possibility to understand better many rel- 
evant phases of matter, such as the celebrated Haldane 
phase of quantum spin chains. First encountered analyti- 
cally at a so-called Affleck-Kennedy-Lieb-Tasaki (AKLT) 
point of the spin-1 antiferromagnetic Heisenberg model 
with a bilinear-biquadratic interaction ' , this phase has 
become a landmark in the study of strongly correlated 
quantum systems in one spatial dimension. 

The existence of a Haldane phase in integer quantum 
spin chains motivated the study of generalizations of the 
spin models with SU(2) symmetry. For instance, gener- 
alized AKLT models with a Valence Bond Solid structure 
have been proposed, which promote SU(2) symmetry to 
some larger Lie groups such as SU(n) 6-8 , Sp(2n) !) and 
SO(n) * together with supersymmetric versions 12 ' 13 . 
Many exciting results also have been found in more gen- 
eral quantum spin chains with such symmetries, includ- 
ing the existence of Haldane-like phases. 

Here we wish to focus on the case of SO(n)-symmetric 
quantum spin chains. In particular, our aim is to further 
contribute to the study and understanding of the dif- 
ferent phases present in the SO(n) bilinear-biquadratic 
quantum spin chain 
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In the above expression L ab (1 < a < b < n) are the 
generators of SO(n) in the n-dimensional vector rep- 
resentation, with the Casimir operator normalized as 
J2 a <b( Lab ) 2 = n - 1, and € [0, 2tt) is the Hamiltonian 
parameter. Despite its simplicity, the model is known to 



have a very complex and rich phase diagram offering a 
wide variety of phases 1 "' 14,15 (see Fig. 1). In fact, this 
phase diagram is not totally understood yet. Even in the 
simple case of n — 3, where the model corresponds to 
the well-known spin-1 bilinear-biquadratic quantum spin 
chain, there is controversy about the presence or not of 
a spin-nematic phase 16-24 near 6 = 5-7r/4 25 . 

In this paper we investigate the phase diagram of 
the above model by focusing on a global property of 
the ground state in the thermodynamic limit, namely, 
the geometric entanglement (GE) . This quantity has 
proven useful in the study of several quantum many- 
body systems and their associated quantum phase tran- 
sitions without the need of order parameters 2 ' -38 . Here, 
we show that for the model in Eq. (1) the GE shows a 
rich behavior as compared to other quantities such as the 
mutual information , fidelity susceptibility : , Renyi en- 
tropies and correlation functions, and is compatible with 
calculations of the fidelity diagram - . More specifi- 
cally, we consider numerically the cases of n = 3, 4 and 5. 
As we shall see, our analysis is able to capture distinctive 
behaviors which we conjecture to be general and valid for 
arbitrary n. For instance, an anomaly in the GE is ob- 
served at the purely bilinear point 9 = 3tt/2, where the 
model is dual to the n 2 -state Potts model and its ground 
state energy and gap can be obtained exactly 49-51 . This 
motivates us to give an intuitive argument in favor of an 
infinite entanglement length and finite correlation length 
at this point' 2 . 

This paper is organized as follows: in Sec. II we explain 
the basic phase diagram of the model for SO(3), and also 
on the known properties of the phase diagram for SO(n). 
In Sec. Ill we comment briefly on the methods we have 
used in our study, including the entanglement measure 
and the numerical technique. Sec. IV contains our re- 
sults for the GE for SO(3), SO(4) and SO(5), together 
with a discussion about these results and the behavior for 
arbitrary n. In this section we also compare our results 
to those obtained by alternative approaches. Finally, in 
Sec. V we briefly summarize our conclusions. For com- 
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pleteness, we include in Appendix A the calculation of the fidelity diagram for the SO (3) case. 




FIG. 1: Phase diagram of the SO(n) bilinear biquadratic chain for (a) SO(3), (b) SO(4) and (c) SO(5). Black and grey dots 
indicate respectively known continuous and discontinuous transitions, whereas white dots indicate other interesting points of 
the system. The white dot close to 9 — corresponds to the MPS point tan 8 = 1/n in all cases (AKLT for n = 3). 



II. BASIC PHASE DIAGRAM 

A. SO(3) 

Let us start by quickly reminding the phase diagram 
of the model for the case of SO (3), see Fig. 1(a). In 
the region 9 G (77r/4,7r/4) the system favors a Haldane 
gapped phase, and at 6 = 0.10247T (tan 6 = 1/3) the 
system corresponds to the AKLT model with an exact 
valence bond ground state 1 ' . At 9 = ir/4 there is a con- 
tinuous phase transition of the Kosterlitz-Thouless (KT) 
type to a critical trimerized phase with k — ±2-7r/3 spin 
quadrupolar correlations, which extends through the re- 
gion 9 € (vr/4, 7r/2). This phase is separated from a gap- 
less ferromagnetic phase in 9 6 (7r/2,57r/4) at a first 
order transition point 9 = ir/2. At 9 = 5ir/4 there 
is another first order transition to a dimerized phase 
that extends through 9 € (57r/4, 7n/4). It has been 
suggested that a spin nematic phase appears at the re- 
gion 9 £ (5tt/4, 1.337r) with a generalized KT transition 
around 9 « 1.337T ' . The existence or not of this spin 
nematic phase is still a matter of discussion . Also, 
many of the properties of the system can be obtained 
exactly at the purely biquadratic point 9 = 3tt/2, where 
it is dual to the 9-state Potts model 49-51 . The quantum 
phase transition at 9 = 77r/4 between the dimerized and 
Haldane phases is of second order. Also, the system is ex- 
actly solvable at the transition points 9 — 7r/4, 7r/2, 57r/4 
and 77r/4. 

B. SO(n) 

The case of general SO(n) is similar, see Fig. l(b,c) 
for n = 4 and 5. In general, the phase that extends 
through the region 9 6 (tan -1 ^-2) 2 > ^ &n ~ 1 r ra -2) ) 



ways contains an MPS point at tan 9 = 1/n. This phase 
is Haldane-like for n = 3, 5 and a non-Haldane spin liq- 
uid for n = 4 53 ' 5 . Field theory analysis seems to indicate 
that there is a KT transition at 9 — tan -1 , 1 „, towards 

(n-2) 

a critical phase 55 . At 9 = ir/2 there is a first order tran- 
sition towards a phase with ferromagnetic order. This 
phase extends up to 9 = tan -1 („1 2 ) 7r ' wnere there is 
a transition to a dimerized phase. Again, at the purely- 
biquadratic point 9 = 3tt/2 the model is dual to the n 2 - 
state Potts model 1 ' 1- ' 1 . We remark that the Heisenberg 
point 9 = does not necessarily reside in the same phase 
as the MPS point. For details see e.g. Ref. 10 ' 14 ' 15 . 



III. METHODS 

A. Geometric entanglement 

Here we choose to explore the details of the above 
phase diagrams for n = 3, 4 and 5 using the geometric 
entanglement per site £ (which we refer to as GE) 26 " 38 . 
Let us briefly remind the basics of this quantity. As- 
sume we are given a quantum state of N parties 
belonging to a Hilbert space n = (g^LjVM, where 
V^l is the Hilbert space of party r. We now consider 
the closest normalized product state of the parties to 
By "closest" we mean the normalized product state 
|$) = I^W) ® |0l 2 l) ® ••• ® \<f>W) that minimizes the 
squared distance |||$) — \^)\\ 2 between |$) and \ty) or, 
equivalently, maximizes the absolute value of their over- 
lap, A max (^) = max$ |($|$)|. This closest product state 
approximation to l^) allows us to quantify its entangle- 
ment via the extensive quantity E(ty) = — log A 2 nax (^'), 
where we have taken the natural logarithm. As de- 
manded, E(^>) is zero if l^) is unentangled. We can also 
define the thermodynamic quantity £ and its finite-size 
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version £pj as 

£ = lim £ N , £ N = N^E^f). (2) 

N— >oo 

The quantity £ in the above equation defines the global 
geometric entanglement per site, or density of global geo- 
metric entanglement. This will be the quantity of inter- 
est in this paper. In the past, this quantity has been 
used in the study of quantum phase transitions 38 , 
quantification of entanglement as a resource for quan- 
tum computation , local state discrimination and, 
more recently, it has also been directly measured 
experimentally . 



B. MPS and iTEBD 

Our results have been obtained by approximating the 
ground state of the model in Eq. (1) by a matrix product 
state (MPS) in the thermodynamic limit 59 ' 60 . As is well 
known, MPSs offer an accurate description of quantum 
states of ID quantum many-body systems, specially for 
gapped systems. Let us be more specific: for a system 
of size N with periodic boundary conditions, these are 
states defined as 

|*)= ]T tv(A^---A^)\z 1 ,...,z N ), (3) 

ii,...,i N 

where A^ lm is a \ x X matrix at site m — 1, . . . , N 
for each i rn — l,...,d m , which labels a local basis of 
the Hilbert space of dimension d m at site m. Since our 
model is an SO(n) spin chain, we have d m = n at all sites. 
Moreover, \ is a refinement parameter for the MPS and 
is called the bond dimension. In practice, MPS offers a 
good variational family of states to approximate quan- 
tum states of many-body systems, specially for ground 
states of Hamiltonians with local interactions in one spa- 
tial dimension. The refinement parameter x control the 
number of variational parameters in the ansatz. There- 
fore, the larger is Xj the more accurate is the variational 
approximation. 

In our case, we choose to use the iTEBD 
algorithm 59 > 60 ( an extension of the MPS method to in- 
finite systems, to obtain approximations to the ground 
state of our quantum spin chain. This method makes use 
of an evolution in imaginary time in order to update the 
different matrices of the MPS at each time step. The 
particular technical details of the method can be found 
in Ref. 59 ' 60 . 

For our purposes, we find that an MPS with inner bond 
dimension ^ up to ~ 120 is sufficient to achieve conver- 
gence in the results' '. Also, the GE is computed from 
the MPS wave function by using similar techniques as 
in 38 , but using product state approximations with a pe- 
riodicity of 6 sites. This helps in avoiding local minima 
in the optimization process to compute the GE. 




FIG. 2: GE of the SO(n) bilinear-biquadratic model, for (a) 
SO(3), (b) SO(4) and (c) SO (5). 



IV. RESULTS 
A. GE for SO(3), SO(4) and SO(5) 

Let us now present our results for the GE of the SO (3) 
model. This is shown in Fig. 2(a) for the whole region 
9 € [0,2tt). Quite remarkably, the GE displays some 
sort of singularity close to all the known transitions in 
the system. In particular, at 9 — 7r/4 it has a cusp, 
which we interpret as an indicator of the KT transition 
between the Haldane and trimerized phases, see Fig. 3(b). 
This is remarkable, since at a KT transition all the two- 
body observables as well as their derivatives are analytic, 
yet the GE is not. A similar cusp in the GE was also 
found at the KT transition of the XXZ spin-1/2 chain 
in Ref. 38 . These cusps in KT transitions are due to a 
sudden change in the closest product state to the ground 
state (see Ref. 38 ). 

Moreover, the GE has a discontinuity close to 9 = n/2, 
which is reminiscent of the first order transition at this 
point. It is also non-analytic at the first order transition 
point 9 = 57r/4. Close to this point it arises from zero 
seemingly according to the power law £ ~ y/6 — 5tt/4, 
see Fig. 3(c). Notice, though, that this continuous behav- 
ior is not incompatible with the presence of a discontinu- 
ous transition at 9 = 57r/4' 4 . We see, thus, that discon- 
tinuous transitions can be detected in the GE as (i) dis- 
continuities (as in 9 = ir/2), or (ii) cusps (as in 9 = 5tt/4). 
Consequently, in order to distinguish if a given cusp in the 
GE corresponds to a discontinuous transition of type-(ii) 
or to a KT transition, one may need to also study e.g. the 
correlation functions of the system or energy derivatives. 
In the first case these quantities would be discontinuous 
across the transition, whereas across the KT transition 
they would be continuous. 

We also observe a peculiar behavior close to the tran- 
sition at 9 = 77r/4, see Fig. 3(f). Specifically, the GE 
begins to drop sharply near 9 w 1.7 Att and begins to 
flatten near 9 ss 1.757T. At 9 w 1.817T the GE sharply 
arises. Indeed, it is possible to estimate the correlation 
length exponent v from the GE close to the transition 
at 9 = 77r/4, since close to a quantum critical point 
£ ~ a\9 c — 9\ v + b . Fitting our results very close to 
this point is compatible with v = 0.5 ± 0.1, in agreement 
with a mean-field theory estimation v = 1/2. However, 
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FIG. 3: (Color online) Different zooms to the GE diagram 
from Fig. 2(a): (a) Around 9 ~ 0.15-7T the GE is smooth (line 
is a guide to the eye); (b) At the KT transition point 9 — n/4 
the GE shows a cusp (line is a guide to the eye); (c) The 
GE goes to zero at the first order transition point 9 — 5n/4, 
seemingly according to the continuous law £ ~ \/(9 — 5n/4) 
(continuous line); (d) At 9 w 1.34n the GE shows a disconti- 
nuity, which eventually vanishes as the MPS bond dimension 
increases (black line is for x — 80, blue line for \ — 100, red 
line for \ — 120, and the position of points has been omitted 
for the sake of clarity); (e) At the integrable point 9 = 3n/2 
the GE shows a cusp (line is a guide to the eye); (f) Around 
the 2nd order transition at 9 — 7n/4 the GE shows a smooth 
change of behavior, compatible with £ ~ a\9 c — $\" with v ~ 1 
(line is a guide to the eye). 



if another fit is done slightly away from 9 = 77r/4 (from 
below), we obtain v ~ 1 instead, as predicted by a con- 
formal field theory approach . This cross-over behavior 
is consistent with the mean-field-like behavior of MPS 
when close to quantum critical points (see Ref. ). 

In addition to the above observations, we see a number 
of other peculiarities in the GE. For instance, the data 
between 9 = 3tt/2 and 9 ~ 1.77T can be very well fitted by 
a straight line. We also observe that near 9 = 0.15-7T it at- 
tains a maximum. A closer look to this maximum shows 
that it is smooth, see Fig. 3(a). The presence of this max- 
imum is intriguing as at this point the ground state pos- 
sesses the highest geometric entanglement in this model, 
albeit we are not aware of this point being special for 
any other reason. Notice also that there is no sign of 
any peculiarity in the GE around the AKLT point at 
9 w 0.10247T. However, we observe two other distinc- 
tive behaviors. First, near 9 w 1.347T the GE seems to 
display a small discontinuity, see Fig. 3(d). As seen in 
the figure, an analysis for large bond dimension shows 
that this discontinuity eventually vanishes as the MPS 
bond dimension increases, and is therefore an artifact of 
the numerical calculation and not caused by the hypoth- 
esized transition towards a spin nematic phase (which 
some studies estimate around 9 rj 1.33tt) » . Second, 
the GE displays a cusp at 9 = 3ir/2, see Fig. 3(e). At 
this point the system is purely biquadratic and dual to 
the nine-state Potts model 49-51 . Such a cusp looks sim- 



ilar to the one in the elusive transition of the deformed 
AKLT model, where the entanglement length diverges 
while the correlation length is finite 38 ' 52 . We will further 
analyze the physics of this point later on. 

In the case of SO(4) and SO(5), the behavior of the GE 
shares many basic features with the previous calculation 
for SO(3), see Fig.2(b,c). Specifically, we observe: (i) a 
peak at the KT transition point 9 = tan -1 ( n _ 2 ) > (ii) a 
discontinuity at the first order transition in 9 = n/2; (iii) 
a non-analyticity at 9 = tan -1 („1 2 ) 11 consistent with 
a discontinuous transition; (iv) a peak at the biquadratic 
point 9 — 37r/2; (v) a linear behavior between 9 = 3tt/2 
and 9 ~ 1.7tt; and (vi) a sudden change of behavior at 
9 = tan -1 (n~2) 2 com P a tible with a continuous transi- 
tion. Notice also that the Heisenberg point 9 = is a 
transition point for SO(4), which is quite different from 
the SO(3) (Haldane) and SO(5) (dimerized) cases. 



B. Discussion and arbitrary-n behavior 

Let us now discuss two aspects of the above results. 
First, notice that all the features observed above for 
SO(4) and SO(5) are also common to our SO(3) cal- 
culations. Thus, we conjecture that these are generic 
properties of the many-body entanglement of the system 
regardless of the value of n. Second, we believe that there 
are good chances that 9 = 3n/2 is a point of infinite en- 
tanglement length and finite correlation length for any n. 

Let us give an intuitive argument in favor of this. 
First, we remind that at this point the ground state sub- 
space is made of two degenerate dimerized states, each 
one of them adiabatically connected (and thus with the 
same long-range properties) to a dimerized state of SU(rt) 
singlets between nearest neighbors, see Ref. ' '"" for a 
proof'* 1 . These dimerized states are known to have infi- 
nite localizable entanglement (and thus infinite entan- 
glement length) if measurements are allowed on pairs 
of nearest-neighbor spins that do not share a singlet. 
Such long-range properties should then also belong to 
the two dimerized ground states at 9 = 3tt/2. Also, since 
this point in the parameter space belongs to a dimer- 
ized gapped phase, translational invariance is broken in 
the ground state of the system towards one of these two 
dimerized states. Therefore, the system should display 
an infinite entanglement length whereas the correlation 
length remains finite because the whole dimerized phase 
is gapped. This, in fact, is very similar to what happens 
in the ground state of the deformed AKLT model' 38 ' 52 and 
two other MPS models (see Ref. 28 ). A quantitative anal- 
ysis of the entanglement length would be needed to ascer- 
tain this intuitive argument. Such an analysis, though, 
is currently beyond our reach. 
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C. Comparison to other approaches 

A comparison between our results and those obtained 
by other methods is in order. It is known that several 
quantities may be singular at the well-known transition 
points of the model, but not around other possible tran- 
sitions. For concreteness, we consider here the case of 
SO(3), and calculations of the mutual information , fi- 
delity susceptibility , fidelity diagram, Renyi entropies 
and correlation functions. Specifically: (i) The mutual 
information was computed in Ref. for finite systems up 
to 14 sites using exact diagonalization, and while it de- 
tects some of the transitions in the system, its behavior is 
analytic around 9 — 7r/4 and 7tt/4, as well as around the 
conjectured nematic phase at 9 = 3ir/2. (ii) In Ref. 40 the 
authors perform a study of the fidelity susceptibility up 
to 12 sites using exact diagonalization, and no anomalous 
behavior is observed at the KT transition at 9 = 7r/4 and 
at the second order transition at 9 = 77r/4. Nothing spe- 
cial is either observed around 9 ss 1.347T and 9 — 3ir/2. 
This fact is interesting since the fidelity susceptibility is, 
essentially, the curvature in the axis direction along a 
diagonal of the fidelity diagram. Its magnitude seems 
then not to be sufficiently strong to identify certain non- 
analiticities in the system that can otherwise be detected 
by looking at the whole picture of this diagram (see sup- 
plementary material), (iii) We have compared our re- 
sults to those obtained by analyzing the Renyi entropies 
of half an infinite chain, which include the von Neumann 
entropy and the single copy entanglement as particular 
cases 63-67 . The Renyi entropies and, in fact, the entan- 
glement spectrum of half an infinite chain 1 ' , are a di- 
rect byproduct of the iTEBD method that we employed 
to approximate the ground state wavefunction. We have 



seen no signatures of anomalous behaviors in these quan- 
tities at 9 = 3tt/2. (iv) Finally, we also computed several 
two-point correlation functions from the obtained MPS 
of the ground state, and saw no clear sign of anomaly 
around the integrable point 9 — 3tt/2. All these observa- 
tions are in contrast with the behavior of the GE. 



V. CONCLUSIONS 

In conclusion, we have investigated the phase diagram 
of the SO(n) bilinear-biquadratic quantum spin chain us- 
ing the GE, and by considering numerically the cases of 
n — 3,4 and 5. We have conjectured that our numerical 
observations are also valid for arbitrary SO(n). Further- 
more, we have seen that the GE provides a remarkably 
rich behavior in the phase diagram as compared to other 
quantities. 

To finish, let us mention that it would be interesting to 
determine which quantities, apart from the GE, can be 
useful to determine the rich phase diagram of the models 
considered here. For instance, a detailed analysis of the 
degeneracies in the entanglement spectrum of the ground 
state for arbitrary SO(n) should be helpful in determin- 
ing the relevant phases' 1,72 . This will be the subject of 
future investigations. 
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This is similar to the Ising chain in parallel field " s . 

75 This follows from Ref. 36 . 

76 These singlets can be interpreted as having SO(n) symme- 
try, but acquire SU(n) symmetry at this point . 



Appendix A: fidelity diagram for SO (3) 

Let us consider the fidelity diagram (FD) of the system 
for the case of SO (3). The aerial view of the (9i,9 2 ) plane 
is represented in Fig. 4 (notice that the FD is symmet- 
ric). As can be seen in the plot, the FD is discontinuous 
at (6*i, #2) = (7r/2, 7r/2) and (57r/4, 57r/4), in accordance 
with the first order transitions present at these points. 
The results indicate that the transition at (tt/2,tt/2) 
seems to be caused by a much stronger discontinuity than 
the one at (5tt/4, 57r/4). Furthermore, a pinch point is 
present at (7r/4,7r/4) corresponding to the KT transition 
in the system. Therefore, these results show that the 
FD is able to capture the presence of this transition, un- 
like other closely-related quantities such as the fidelity 
susceptibility. The presence of another pinch point at 
(7tt/4, 7tt/4) is compatible with the second order quan- 
tum phase transition here. 

Importantly, we see a number of other features in 
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the FD. First, we observe a tiny pinch point around 
« (1.347T, 1.347r). This would be compatible with a tran- 
sition towards a spin nematic phase at this point, yet we 
stress that we can not conclude on the existence of such a 
transition based on these results only. Second, we observe 
a line of low fidelity between the points PI = (7r/4, 57r/4) 
and P2 = (tt/2, 3ir/2) in the lower half of the plane (and 
its symmetric counterpart in the upper half plane) . There 
are two remarkable facts about this line. First, notice 
that it goes between points that are known to be spe- 
cial in the model: at 9 = 7r/4 there is a KT transition, 
at 9 = tt/2 and 9 = 5tt/A there are first order tran- 



sitions, and at 9 = 3ir/2 the model is exactly solvable. 
Intriguingly, this last one is the parameter value at which 
the GE shows a cusp, see Fig. 3(e). Thus, it looks like 
9 = 37r/2 is also special from the point of view of the 
FD. Second, the line of low fidelity is slightly interrupted 
around P3 ~ (1.347T, 0.377r). This is another indication 
that 9 ~ 1.34 is special for the FD (notice also that 
9 = 0.377T is not special: it is simply the y-coordinate 
over the line between PI and P2 when the x-coordinate 
is ~ 1.34). This behavior is also in accordance with our 
previous observations using the GE. 



